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Abstract. We discuss the static and dynamic multicritical behavior of three-dimensional systems of 
© 0{n±) symmetry as it is explained by the field theoretical renormalization group method. 
Whereas the static renormalization group functions are currently know within high order expansions, 
we show that an account of two loop contributions refined by an appropriate resummation technique 
gives an accurate quantitative description of the multicritical behavior. One of the essential features 
of the static multicritical behavior obtained already in two loop order for the interesting case of an 
antiferromagnet in a magnetic field (ny = 1, = 2) are the stability of the biconical fixed point and 
the neighborhood of the stability border lines to the other fixed points leading to very small transient 
exponents. We further pursue an analysis of dynamical multicritical behavior choosing different forms 
of critical dynamics and calculating asymptotic and effective dynamical exponents within the minimal 
subtraction scheme. 
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INTRODUCTION 

Beneath the milestone contributions of N.N. Bogolyubov that shaped modern theoretical physics 
one definitely should mention his and D.N. Shirkov work on the renormalization group (RG)[1]. 
Three papers on RG written in the mid-50-ies by three different groups [2] addressed quantum 
electrodynamics problems, however very soon their importance has been realized in - on the first 
sight - very different field of phase transitions and critical phenomena. It is generally recognized 
by now that the success in conceptual understanding and quantitative description of behavior in the 
vicinity of critical points in different condensed matter systems is due to the effective application 
of the RG ideas originating from the above papers [3]. It is our pleasure to contribute to these 
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FIGURE 1. Typical phase diagrams of anisotropic antifenomagnets in a uniform parallel external mag- 
netic field H. Types of ordering are schematically shown by arrows, a: the bicritical point. Three phases 

- an antiferromagnetic phase, a spin flop phase and the paramagnetic phase are in coexistence. The phase 
transition lines to the paramagnetic phase are second order transition lines, whereas the transition line 
between the spin flop and the antiferromagnetic phase is of first order, b: the tetracritical point Four phases 

- an antiferromagnetic phase, a spin flop phase, an intermediate or mixed phase and the paramagnetic phase 

- are in coexistence. All transition lines are of the second order in this case. Also indicated is the dynamical 
universality class of the transition from the paramagnetic to the corresponding ordered phase according to 
the classification of Hohenberg and Halperin jlSll for the three component antiferromagnet. 



Proceeding^ by a short review of recent work done by application of the field theoretical RG 
approach to analysis of multicritical phenomena. 

Multicritical points appear on phase diagrams of various systems that contain several phase 
transitions lines. In the vicinity of the meeting points of such lines the multicritical behavior is 
observed, which is characterized by competition of different types of ordering. Prominent examples 
are given by the antiferromagnets in an external magnetic field like GdALOs, MnF2, MnCl24D20, 
Mn2AS4 (A=Si or Ge) ft. Other examples are given by the layered cuprate antiferromagnets like 
(Ca,La)i4Cu2404i. Schematic phase diagrams of such systems are shown in FigHJin aH-T plane. 
There, multicritical points of two different types are manifested. At a bicritical point (Fig. [1^) 
three phases are in coexistence, whereas four phases coexist in the tetracritical point (Fig. [D)). 
On a more general level, the multicritical behavior is inherent to a critical system when some 
"nonordering" field is applied. Such a field (beside the magnetic field H this may be pressure, 
stress etc.) may alter non-universal parameters of the system and lead to appearance of the lines 
of phase transition points. Besides the above example that concern the shift of the Neel point of 
anisotropic antiferromagnets by a uniform magnetic filed, other examples of multicritical behavior 
are observed at a shift of the Curie points under applied pressure or depression of the A point in 
^He at dilution by ^He fl. 

A field theoretic description of multicritical behavior starts with a static effective Hamiltonian 
for an n-component field 4> — (^||, of © 0{n^) symmetry (n|| -\-n^ = n). An account of 
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the interaction between the two order parameters ^|| and ^± leads to different types of multicritical 
behavior connected with the stable fixed point (FP) found in the RG treatment 10, 0, fl 0, 



111 uM- In particular, the bicritical point (Fig. [T^) has been connected with the stability of the 
isotropic Heisenberg fixed point of 0{n\\ + n^) symmetry, whereas the tetracritical point (Fig. [l]b) 
corresponds to a FP of 0{n\^\) © 0{nx) symmetry, which might be either the so called biconical 
FP or the decoupling FP. In the last FP the parallel and the perpendicular components of the order 
parameter are asymptotically decoupled. If no FP is reached the multicritical point might be of first 
order, i.e. a triple point. 

Quite recently the possible types of phase diagrams in the H — T plane of three dimensional 



uniaxial anisotropic antiferromagnets have been studied by Monte Carlo simulations [|l2h . For 



n|| = 1 and n \ = 2 a phase diagram with a bicritical point has been found in agreement with earlier 
simulations [ll3|], but contrary to the results of RG theory in higher loop orders duj]. 

The dynamics of antiferromagnets in a magnetic field is quite complicated. To account for 
conservation laws present in such systems, the dynamical equations of motion should contain 
coupling terms between the two order parameters (the components of the staggered magnetization 
parallel and perpendicular to the magnetic field, ^|| and and conserved densities (e.g. the 
parallel component of the magnetization or the energy density). First formulation of the equations 
of motion at multicritical points has been done in Ref. ||l7l]. The simplest form of dynamics 
assumes the relaxational behavior for the two order parameters ^ii and (the so-called model A) 
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18ll . Dynamical multicritical behavior within the one-loop approximation has been considered 
in [|l7|l on the basis of the static one loop results [|9|]. A further step to the complete model is to 
include the diffusive dynamics of the slow conserved density_leading to a model C like extension. 
This model has been studied in one loop order in Refs. HviSH] taking into account only a part 
of dynamical two loop order terms and one loop statics. In order to get more insight in the dynamics 
in the vicinity of multicritical points, recently we have reconsidered the above dynamical models 
within the two loop approximation [fl^llll]. 

In what follows below we briefly summarize an outcome of an RG analysis of multicritical 
behavior paying special attention to an impact of the non-universal contributions to an asymptotic 
behavior. In particular, we will show that an account of two loop part of the RG expansions refined 
by an appropriate resummation technique gives an accurate quantitative description of the static 
multicritical behavior. Furthermore, we pursue an analysis of dynamical multicritical behavior 
choosing different forms of critical dynamics and calculating asymptotic and effective dynamical 
exponents. 

RG FLOWS AND STATIC MULTICRITICAL BEHAVIOR 

The generalized static 0{n\\) © C)(nx) -symmetrical effective Hamiltonian that results from the 
decomposition of the n-component order parameter field into two mutually interacting fields ^|| 
and of different irreducible representations of dimensions n|| and n^, n — n\\^-\- nj_, reads: 

f fl- - 1"^- - 1- - 
^ J ^''^ j 2^"^^^^ • ^±0 + 2 L ^'^^0 ■ ^4±o + 2''||^||o ■ ^||o + 2 E ^'^l|o ' '^'^||o 

+ ^(^±o-^i.o) +^(^||o-^||o) + ^ (^±0 • ^j.o) (^||0 • ^||o) ^ ■ (1) 



Here, {u±jUx,u^^} = {u} and r±, f|| are couplings and masses, correspondingly, index refers 
to the bare quantities, and central dots stand for scalar products. The decomposition in parallel 
and perpendicular order parameter components allows to describe the multicritical behavior at the 
meeting point of two critical lines: (i) the line where r± becomes zero and the n^-dimensional 
components ^j^o are the order parameter, and (ii) the line where becomes zero and the order 
parameter is ^| q- At the meeting point both quadratic terms become zero and both components 
of have to be taken into account. As has been predicted already by the one-loop RG analysis 
iSImI, an effective Hamiltonian ([T]) describes three different types of multicritical behavior that are 
governed by three different FPs: (i) the isotropic n component Heisenberg FP, called below J^{n), 
all fourth order couplings are equal in this FP, (ii) the decoupling FP point ^, which consists of a 
combination of the FPs J^{n±) and J^(n||) of two decoupled systems and (iii) the biconical FP, 
with nontrivial nonzero couplings. As it was revealed by subsequent calculations CiEll] the 
FP picture does not change qualitatively with an account of higher orders of the perturbation theory. 
However, the one-loop results attain essential quantitative changes that lead to drastic modification 
of the type of a phase diagram. Typical example may be given by the behavior of the j8-functions, 
that describe flow of the fourth order couplings {u} under renormalization. The above functions. 



calculated in the two-loop approximation with the minimal subtraction RG scheme read [|l4|l 



o (n^ + 8) 2 "II 2 (3n± + 14) ^ 5m n,, 

= -£u^ + ^^^ul + -^ul-^^ '-ui-^u^ul-^ui, (2) 

„ , (nx + 2) , («||+2) 2 2 (nx + n|| + 16) 



(^± + 2) 2 h+2) 2 5(n^ + 2) 2 5(n||+2) ^ ... 

- UxU± Mx"|| "±"x MxM|[, (3) 

R I ("11 + ^^ ,2 I 7 (3^11 + 14) 3 5n_L 2 3 ... 

P«li = -ew|iH ^ — "|| + -^"x Y2 "II ~ ^^"ll"x - -^"x- (4) 

Here, {u±,Ux,u\\} = {u} are renormalized couplings and the space dimension d enters the /3- 
functions via parameter £ = 4 — d. With the /3 -functions at hand, one can analyze the flow equations 
of the fourth-order couplings {u}: 

^^ = AJM), (5) 

with a =±, II , X and the flow parameter £, and find the FPs {u*} of these equations as the solutions 
of the system of equations 

AJK}) = o. (6) 

Which of these FPs is the stable one depends on the number of components n± and n|| and the 
dimension d of space. The scaling properties depend on the symmetry of stable FP. 

There are two alternative ways to look for the solutions of the FP equations Q and, subsequently, 
for the scaling properties of the system. In one approach, the £-expansion, the solutions are 
obtained as series in e and then evaluated at the value of interest (at e = 1 for d = 3 theories). 
Alternatively, one may solve a system of non-linear equations directly at the dimensionality of 
space of interest (e.g. at £ = 1) [22] and obtain the FP coordinates numerically. The RG expansions 
being divergent [23], the special resummation techniques are used to get convergent results [24]. 
As we have discussed already above, depending on the values of n^, and d, the multicritical 
behavior is governed by one of the three non-trivial FPs: {u*j_ = u*x = u*^}, ^ {m^ 7^ m* 7^ 
and '3) \u\ 7^ 0,m* = 0,Mj| 7^ 0}. In Fig. |2] we show how the stability of these FPs change with 
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FIGURE 2. Regions of FPs stability in the in the n|| — -plane, d = 3. The lines separate regions where 
Heisenberg FP J^, biconical FP and decoupling FP ^ are stable (from left to right). Shown are the 
J^^-stability borderUnes (dashed lines) and =^^-stability borderlines (soUd lines), in one loop order (thin 
lines) and resummed two loop order (thick lines). The dots indicate low integer values for order parameter 
components. 



n±, for d = 3. There, we compare the first order e-expansion results [jSi |91] with the two-loop 
results il4 \ obtained within the fixed d = 3 technique 12211 . The two-loop results were obtained 
applying Pad-Borel resummation technique to functions (I2l)-(|4l) [251]. One sees that the borderlines 
of the FPs stability are drastically shifted to smaller values of OP components. Thus in the case 
n| = 1 and = 2 FP ^ (connected with tetracriticality) is stable in two loop order contrary to 
the one loop calculations where the FP J^f (connected with bicriticality) is stable. The resummed 
higher orders of the perturbation theory do not change this result and do not lead to essential 
changes in the critical exponents either [ jlll ]. 

As usually, the asymptotic values of the critical exponents are defined by the stable FP values of 
the corresponding RG ^-functions, which we do not expose here. Note that in general, there are 
distinct exponents 77||, r}± governing spacial decay of the order parameter correlations in directions 
parallel and perpendicular to the anisotropy axis. As a consequence, there is a pair of y-exponents, 
yj|, 7x that govern corresponding isothermal magnetic susceptibilities. However, the above RG 
procedure assumes that the multicritical system is described by a single diverging length scale and 
therefore by one correlation length ^ and one corresponding critical exponent v. This does not 
hold for decoupled systems where two length scales are present and the usual scaling laws with 
one length scale break down [9]. We give typical numerical values of the exponents in Table [U 

Whereas the asymptotic critical exponent values are determined strictly at the FP and correspond 
to the scaling behavior at the multicritical point, of special interest are the effective critical 
exponents which are observed in the vicinity of the multicritical point. These are the effective 
exponents that often are observed experimentally and are measured in MC simulations. In the 



TABLE 1. Critical exponents of the 0(1) © 0(2) model obtained in different approx- 
imations. 1 141: resummation of the two-loop RG series at fixed d = 3; |@]: first order 



e-expansion, Bill 12611 : resummed fifth order e-expansion. Numbers, shown in italic were 
obtained via familiar scaling relations. 
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RG framework, one may estimate the effective exponents from the values of corresponding RG 
^-functions calculated along the RG flow and relate the flow parameter i to the distance to the 
multicritical point. In Fig. [3] we show the resummed [^25^ RG flow of Eqs. ([5]) for different initial 
conditions The unstable FPs are shown as filled spheres, the stable biconical FP as filled cube. 
Let us note that the neighborhood of the stability border lines to the other FPs leads to very small 
transient exponents. Therefore, the stable FP is not reached for the value of the flow parameter 
chosen in Fig. [3] (there, the flow parameter has been changed in the interval —40 < \ni < 0). 




FIGURE 3. Resummed RG flow of Eqs. ^ for different initial conditions at d = 3, = I, n± = 2. The 
unstable FPs are shown as filled spheres, the stable biconical FP as filled cube. The FPs points are connected 
by separatrices defining the surface which encloses the attraction region. 

Defining the effective exponents as explained above, one can evaluate their numerical values 
along the RG flows of Fig. [3] and in this way predict possible outcome of measuring the scaling 
properties of different observables at multicritical point. As two typical examples, we show in 
Fig. HI the change of the values of isothermal susceptibility effective exponents yj|, Y±^ and of 
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a. b. 

FIGURE 4. Effective exponents of different observables in the vicinity of a multicritical point for the the 
flows of Fig. |3] a: isothermal magnetic susceptibility (solid curves: ]^|, dashed curves: y^)', b: correlation 
length. 



the correlation length critical exponent v as the multicritical point is being approached, the limit 
T Tc corresponds to the limit £ ^ 0. 

Before passing to discussion of peculiarities of dynamic multicritical behavior, let us note a 
particular feature of the © 0{n^) model that becomes evident from the above analysis of the 
statics. As the stability analysis shows, for the physically interesting case d = 3,n\\ = l,n^ = 2 the 
asymptotic behavior is governed by the biconical FP Therefore, the tetracritical point is realized 
(c.f. Fig. [T^). However, depending on the particular microscopic non-universal characteristics of a 
given system, one may expect a variety of different scenarios for multicritical behavior, including 
the triple point (that corresponds to the run away solutions of the RG flow equations, c.f. Fig. [3]) 
and bicritical point (when for certain initial condition the Heisenberg FP is reached). 

DYNAMICS IN THE VICINITY OF MULTICRITICAL POINTS 

Sketched above particular features of static multicritical behavior are further manifested if the crit- 
ical dynamics is addressed. Below, we briefly analyze three different from of dynamical behavior 
in the vicinity of multicritical points. 



Relaxational dynamics (model A) 

Let us start from the simplest dynamical model, model A, when one assumes relaxational 
behavior for the two order parameters ^|| and ^j^. This model has been studied in the one-loop 
approximation in [17], the two-loop results have been obtained in [14]. The model A type Langevin 
equations of motion describe two order parameters that relax to equilibrium with the relaxation 
rates (kinetic coefficients) T±_ and r||: 

- + (7) 



dt 



-r 



5^ 



(8) 



Here, J^f is the static effective Hamiltonian ©, index refers to bare (unrenormalized) quantities 
and the stochastic forces Otj, , , Oa. fulfill Einstein relations 



(0^«(x,O<(V,O) = 2t^d{x-x')dit-t')5''P 
= 2^d{x-x')d{t-t')d'J , 



(9) 
(10) 



with indices Cf, j8 = 1 , . . . , and /, 7" = 1 , . . . , n|| corresponding to the two subspaces. 

Application of the RG procedure to study dynamical multicritical behavior relies on the Bausch- 
Jans sen- Wagner approach fi27i1 . where the appropriate Lagrangian of the model is studied and 
dynamic vertex functions are calculated in perturbation theory and renormalized. In such a tech- 
nique, essential simplification of calculations is achieved due to the possibility to single out a static 
part of every dynamic vertex function [|28l. 12911 . Renormalization of the kinetic coefficients gives 
rise to appropriate j8 -functions. Here, we reveal the two-loop j8 -function for the time-scale ratio 
V = r |/r^ between the renormalized kinetic coefficients Fii and F^. The function reads 1II6II : 



j8v = :^<| [(n||+2)4-(n^+2)Mi](61n^-i; 



■nil u. 



V 2 + v v(2 + v) 



4vl„?(i±^ + 21„M-l 



l+2v 



\+2v 



(11) 



As we have noted in the preceding section discussing the static critical behavior, a non universal 
effective critical behavior may be observed if the values of the static couplings and the time scale 
ratio are not in a FP but rather are described by the flow equations. For v the flow equation reads 



„dv 



(12) 



Below we will show some results about non-universal dynamic multicritical behavior obtained 
with two-loop accuracy. The numerical results for the static part of the RG function were obtained 
by means of the resummation technique n25\\ . whereas no resummation has been applied to the 
dynamic functions 1II6I1 . 

One of the quantities of interest that characterize dynamic critical phenomena is the autocorrela- 
tion time T. It is known to diverge as the critical point Tc is approached, the divergency is described 
by the power law: 

T-lr-Tcp^'^ (13) 



with the universal correlation length and dynamic critical exponents v and z, correspondingly. In 
the multicritical phenomena we consider, one distinguishes two dynamical critical exponents, z\\ 
and z±, that govern the power law increase of the autocorrelation time for the order parameters 
^11 and correspondingly. In asymptotics they are defined by the stable FP values of the 
corresponding RG functions. At the strong scaling FP there is only one dynamic time scale and 
the two exponents are equal whereas at the weak scaling FP they are different and define for each 
component, parallel and perpendicular, the time scale. As it follows from our calculations [16] and 
as one may see from the Fig. [51 the region of stability of the biconical FP ^ (physically important 




FIGURE 5. Regions of the different types of dynamic scaling behavior, e = 4 — d = I. The rest of 
notations aie as in Fig.[2l 




a. b. 

FIGURE 6. Model A multicritical dynamics. Effective dynamical exponent in for different RG flows in 
the vicinity of a multicritical point at d = 3, n^^ = I, n± = 2. The labeling of the flows corresponds to Fig. 
[3l The exponents for the perpendicular (dashed curves) and parallel (solid curves) components of the OP 
differ in the non asymptotic region. 



case J = 3, n|| = 1, = 2 belongs to this region) is characterized by the strong scaling dynamics: 
the time relaxation of both order parameters, ^|| and ^± is governed by the same exponent. In Fig. 
[6] we show an evolution of this exponent Zeff to its asymptotic value z = 2.05 when the time- scale 
ratio V is set to its FP value and the static couplings u change along the RG flows of Fig. [31 Since 
the exponents have not reached their (equal) asymptotic values differences between the parallel 
and perpendicular components of the OP remain. 



Conservation of magnetization (model C) 



A step towards making the description of dynamic phenomena in the vicinity of a multicritical 
point more realistic is to take into account possible couplings between the order parameters 



and conserved densities, that is to consider the model C dynamics 11151 11811 . In the problem 
under consideration, there are two types of conserved densities: one is magnetization-like (more 
precisely, it is the parallel component of the magnetization), another is the energy density. We will 
not consider this second density here, as far as up to the two-loop order the specific heat critical 
exponent a is negative for the case d = 3,n\\ = l,n^ = 2 which is of most interest here. Therefore, 
a coupling to the energy density is irrelevant in the RG sense - it vanishes at the FP [29]. An 
account of both the order parameter and the (conserved) scalar density is achieved by an extension 
of the static functional ([I])- Now, the corresponding model C static functional reads: 

Here, the first term in the right hand side is given by Eq. ([1]), the density mo = mo(x) is a scalar 
quantity, /i is a field conjugated to mo, 7_l and fj| are asymmetric static couplings between the 
corresponding order parameters and the conserved density. 

In their turn, the relaxational equations of motion ([T]),® are now extended by including a 
diffusion equation for the scalar density: 

dt 5(^)||o "^'i 

= AV2— + 0,„. (17) 

ot oniQ 

Here, the static functional J^^^^ is given by (fT4l) . A is a kinetic coefficient of diffusive type for 
the scalar density, the rest of notations is as in ([T]),®. The stochastic forces 0^^, Otp^^ satisfy the 
Einstein relations (|9l), (fTOl) . with an additional Einstein relation for the new stochastic force 6m' 

{dm{x,t) dm{x\t'))=-2Xw^5{x-x')5{t-t'). (18) 

The renormalization of the above introduced asymmetric couplings y^, fj and kinetic coefficient 
A leads to new RG functions. In particular the RG flow of the time scale ratios 

r± Til 
= T ' ^11 = T ^^^^ 

is now governed by the appropriate functions jS^^ and ^ correspondingly. Note that defined for 
model A time scale ratio v is equally well defined in terms of (fT9l) : 



Fii wii 

v=-^ = -^. (20) 

r_L W_L 

Therefore, the dynamical FP equations: 

/3,Jwl,^,v*) =/3,,„(wl,w|,v*) =/3,(wl,^,v*) =0 (21) 



n_L = 2.0 




a. b. 

FIGURE 7. Model C effective dynamical multicritical behavior at J = 3, = 1, «x = 2. a: dynamical 
RG flow for different initial conditions numbered from 1 to 4. b: effective dynamical exponents Z||, zj^, and 
Zm calculated along the RG flows of Fig. |7^, as indicated by the numbers. 



are now not independent: one of these equations can be eliminated by the relation (|20l) . 

Equations of motion (fT5])-(fT7]) describe time evolution of three different observables. Each of 
them has its own autocorrelation time which, as the multicritical point is reached, may be governed 
by an independent dynamical critical exponent. In addition to the two exponents defined in the 
former subsection, Z|| and z±, the dynamical critical exponent Zm for the scalar density is to be 
considered. Similar, as in the model A case, these three exponents may coincide, in the strong 
scaling dynamical FP or they may differ, in the weak scaling dynamical FP. Complete stability 
analysis of the model C RG equations in two-loop approximation is given in Ref. nlw . In particular, 
it is shown that for the case d = 3, = I, = 2, where the static FP is the biconical FP the 
strong scaling dynamical FP is stable. Physically this means that in asymptotics the multicritical 
dynamics is characterized by one time scale, and three dynamical exponents coincide. In particular, 
their asymptotical value was found to be = = = 2. 18 u4\ . However, as it was revealed in 
the former sections, the effective multicritical behavior is much reacher. In particular, in Fig.|7h we 
show the RG flows calculated for different dynamical initial conditions when the static couplings 
are chosen to be fixed at their biconical FP values. The stable dynamical (strong scaling) FP lies 
outside the region shown. Also shown is the surface v = w\\/w± to which the flow is restricted 
by the condition (|20l) . The RG flows of Fig. |7^ give rise to difference in the effective dynamical 
critical exponents, as shown in Fig.|7]b. The insert of the figure shows that even for flow parameters 
as small as \ni = —2000 the effective exponent z± has not reached its asymptotic value 2. 18. 



The complete dynamic model (model G) 

We now restrict ourselves to the case of n|| = 1, = 2 and include mode coupling terms, which 
correspond to Larmor terms describing the precession of the alternating magnetization and the 
magnetization around each other. They are well known from the isotropic antiferromagnet without 



an external field OOh . Then within an external magnetic field the corresponding equations read 

^ = iv^^+^.»<',5.„^+e„. (24, 

Now a and /3 indicates the planar components x,y and the Levi-Civita tensor £"^~ with the third 
index fixed to z has been introduced. The parallel component of the OP is its z-component. This 
component remains just relaxing, whereas the planar components of the OP are coupled to the 
z-component of the magnetization by the precession terms. 

A new feature arises because of the simultaneous presence of the mode coupling g and the 
asymmetric static couplings Y± and in J^^'^^ ([14)) . The perpendicular relaxation coefficient f ^ 
has to be considered a complex quantity where the imaginary part constitute a precession term 
(second term on the right hand side of (l22l) . Even if in the background such terms are absent they 
are produced by the renormalization procedure. 

The stochastic forces 0^^, O^j,^^ and fulfill Einstein relations 

(d^^ixj) ePjx'j')) = 2f'^5(x-y)5(?-05«^ , (25) 
0<|,||(^',O) = 2^5ix-x')5{t-t') , (26) 
{e„,{x,t) Omix'j')) = -2lw^5{x-x')5{t-t') . (27) 



This model has been solved in one loop order in [jlTI] using the one loop results of statics. As is 



has been already seen for the simpler dynamic models models changes are expected in two loop 
order both by the statics as well as by the dynamic terms especially of model C type. We have 
calculated the complete field theoretic functions in two loop order [|3l|] necessary to calculate the 
critical (effective) dynamical exponents. Independent whether the Heisenberg or biconical is the 
stable static FP a first inspection of the flow of the dynamical parameters shows the following: (i) 
The imaginary part of the perpendicular relaxation rate renormalizes to zero, (ii) the times scale 
ratios v (|20l) . w\\ (fT9l) approach zero and w±_ increases to oo. Irrespective of the kind of the stable 
dynamic FP - wether it is a strong scaling FP with very small but finite or a weak scaling FP with 
zero values for v and W|| - the physical observable features of the magnetic transport coefficient 
are effective ones. The range of effective values for the dynamic exponents corresponding to the 
relaxation of the perpendicular and parallel alternating magnetization and the magnetization are 
starting around its Van Hove values ~ Z|| ~ z;^ ~ 2 in the background and approach for the 
biconical FP deep in the asymptotic regime 

zx~1.6 ^11 ~ 2 z;(,~1.6. (28) 

The main prediction according to this result would be that the perpendicular and the parallel 
component of the OP would scale differently in this region. 

The importance of this magnetic system lies in the physical accessibility of the OP, contrary to 
superfluid ^He or superfluid mixture of ^He and ^He whose dynamics is described by model F 



Il32ll . Here all quantities are in principle measurable quantities. Thus the prediction of the different 
dynamic scaling of the OP components can be tested. 



CONCLUSIONS AND OUTLOOK 

By this review we wanted to summarize recent progress achieved in theoretical description of 
the multicritical phenomena. Whereas traditionally RG techniques address critical points in their 
different realizations, the description of multicritical phenomena is possible both on quantitative 
and accurate qualitative levels. Moreover, the problem appears to be tractable analytically even if 
the complicated forms of multicritical dynamics are confronted. As is revealed by the theoretical 
analysis, a particular feature of static and dynamic behavior inherent to multicritical points is the 
multitude of fixed points that describe the RG flow. In its turn, this gives rise to rich effective 
behavior that may be characterized by different types of multicritical points. A natural continuation 
of performed studies would be to analyze cumulative effects caused on the multicritical behavior 
by symmetry breaking factors of different forms (single-ion anisotropics, disorder, frustrations) 
that might be present in a system. 
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